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Abstract 



< 

. , The classical model of q-damped oscillator is introduced and solved in 

_^ ■ terms of Jackson q-exponential function for three different cases, under- 

damped, over-damped and the critical one. It is shown that in all three 
cases solution is oscillating in time but is unbounded and non-periodic. 
By q-periodic function modulation, the self-similar microstructure of the 
solution for small time intervals is derived. In the critical case with de- 
generate roots, the second linearly independent solution is obtained as a 
^ ' limiting case of two infinitesimally close roots. It appears as standard 

derivative of q-exponential and is rewritten in terms of the q-logarithmic 
function. We extend our result by constructing n linearly independent set 
of solutions to a generic constant coefficient q-difference equation degree 
N with n degenerate roots. 



1 Introduction 

The damped harmonic oscillator as simplest classical model of motion with dis- 
V~j ■ sipation, corresponds to friction force proportional to velocity of motion. It 

r% \ appears in many physical problems from quantum theory to inflating universe 

models. The quantum damped oscillator as one of the simplest quantum sys- 
tem displaying the energy dissipation, has been studied to understand dissipa- 
tion in quantum theory pQ. Most popular models are, the Bateman-Feshbach- 
Tikochinsky oscillator as a closed system with two degrees of freedom and the 
Caldirola-Kanai oscillator as an open system with one degree of freedom and 
time dependent mass. Due to complicated character of the friction force, sev- 
eral modifications of the damping term were proposed as fractional derivative, 
time-delay or finite difference derivative etc. The goal of the present paper is 
to study classical q-extended damped harmonic oscillator, where standard time 
derivatives are replaced by Jackson q-derivatives [?. The q-extension of har- 
monic oscillator and its solution in form of basic trigonometric functions [I] was 



considered before in [5] . The q-extended heat equation and corresponding Burg- 
ers equation with q-shock solitons were studied in [5]. In the limit q — > 1, the 
q-deformed model reduces to the standard damped oscillator model. We con- 
struct solution in terms of Jackson q-exponential function and find q-periodic 
modulation of the solution with self-similar properties. Special attention is 
paid for degenerate roots case. And results are generalized for arbitrary order 
constant coefficient q-ordinary difference equation. This gives background for 
further possible quantization of corresponding model. 

2 Damped Oscillator 

In reality a spring never oscillates forever, since frictional forces will diminish 
the amplitude of oscillation until the rest. In many situations the frictional force 
is proportional to the velocity of the mass as follows /,- = — jv, where 7 > 
is the damping constant. Therefore, by adding this frictional force we have the 
following equation for a spring 



d 2 x dx 



in 



dt 2 ' dt 



7— + kx = 0. (1) 



Solution of this equation in the form x(t) = e xt leads to the characteristic 
equation 

m\ 2 + 7A + k = (2) 

with two roots 

—7 + \/^f 2 — Amk —7 — \Z"f 2 — 4mk 

Ai = , A2 = • 

2m 2m 

Then according to value of damping constant we have three cases : 

i- Under-damping Case: When 7 2 < 4mk, which means that friction is 

sufficiently weak, we have two complex conjugate roots 

Al < 2 = -i ±?W ' (3) 



where ui = \J — — -^j ■ Then the general solution of ([T]) is 

x(t) — e~^* (AcosLot + Bs'mujt) . (4) 

If 7 = 0, there is no decay and the spring oscillates forever. If 7 is big, the 
amplitude of oscillations decays very fast (the exponential decay). 
ii- Over-damping Case: When j 2 > Amk, which means that friction is suf- 
ficiently strong, both roots are real, this why the solution decays exponentially 



x(t) = Ae ^ l + Be ^ . (5) 



This case is called as over-damping because there is no any oscillation. 
iii- Critical Case: For 7 2 = 4mk, we have two degenerate roots 



7 
2m' 



Ai — A2 

then the general solution is 

x(t) = Ae-& t + Bte-&*. (6) 

3 q-Harmonic Oscillator 

Here we introduce the q-Harmonic oscillator. Equation of q-dcformed classical 
harmonic oscillator is 

D q x(t)+uj 2 x(t)=0, (7) 

where the q-derivative is definite as |2], 

_ x(qt) - x(t) 
DqX{t) ~ (q-l)t ■ (8) 

Using the power series method (or the q-exponential form x(t) = e q (Xt)), we 
find the general solution of g-Harmonic Oscillator in the following form [3] 

x(t) = A(t)cos q ujt + B(t)am q ujt, (9) 

where 

D q A(t) = D q B(t) = 0, 

means A(t),B(t) in general are g-periodic functions, and particularly could be 
arbitrary constants. Here the Jackson q-exponential function is definite as 

00 ±n 

n— 1 L Jy 

and 

e g (ii) = cos 9 t + isin q t, (11) 

where [n] q = l + g + ... + g" _1 . For g > 1 this function is entire analytic function, 
so we restrict consideration by this case only. 

In Figure[T]we plot particular cos g t solution of g-deformed classical harmonic 
oscillator. In contrast to standard sin t and cos t functions, sin g t and cos 9 t func- 
tions [1], are not bounded and also have no periodicity. In Figure[5]we plot mod- 
ulation of the same solution with q-periodic function A(i) — sin I j^l m t J cos g t, 
which gives micro oscillations to the solution. 
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Figure 1: q-Harmonic oscillator solution cos q t 
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Figure 2: q-Harmonic oscillator solution sin I ^- mil cos q t 



4 q-Damped Harmonic Oscillator 

We define equation for (/-analogue of damped oscillator in the form 



where 



D 2 x(t) + TD q x{t) + oo 2 x{t) = 0, 



V m m 



By substituting x(t) = e q {Xt) into equation (|12|) . we obtain 

e q (Xt) [A 2 +TA + o; 2 ] =0. 



(12) 



(13) 



For q > 1, e q (\t) is an entire function defined for any t, this why it has an 
infinite set of zeros (no poles) . Then, we can choose 

X 2 +TX + LU 2 = 0. 

The roots of this characteristic equation are 

r 



Ai. 



± 



£2 

T 



W z 



4.1 Under-Damping Case 

For F 2 < 4w 2 , we have two complex conjugate roots 



where 



Ai = --+ifl, 



n 



\o = — — — iCl, 



£2 

T' 



Then the general solution of equation (fT2"j) is 

r 



x(t) = Ae q 



iVL)t 



Be n 



m t 



(14) 



In Figure [3] and Figure S] we plot particular solutions with constant (^4 
B = 1) and with q-Periodic modulation, respectively. 
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Figure 3: Under-damping case A = B = 1 



4.2 Over-Damping Case 

For F 2 > 4a; 2 , we have two distinct real roots Ai % and solution is 



x{t) = A{t)e q 



r 

2 



?--' 



B{t)e q 



-\-r-i--v 



(15) 
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Figure 4: Under- damping case with g-periodic function 



where A(t),B(t) are g-periodic functions (or could be arbitrary constants). 

In Figure [5] and Figure [6] we plot particular solutions with constant (A 
B = 1) and with g-Periodic modulation, respectively. 




Figure 5: Over-damping case A = B = 1 



4.3 Critical Case 

For T 2 = Au 2 



The hrst obvious solution 



we have degenerate roots Ai,2 = — §• 
Lot). However if we try the second linearly independent solution in the 



is e q 

usual form te q {~ujt), it doesn't work. This why we follow the next method: 
We suppose that the system is very close to the critical case so that ^ = uj- 
where e <C 1 . Then the roots of characteristic equation are 



Ai 



A, 



(16) 
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Figure 6: Over-damping case with q-periodic function 



and the solution is 

x(t) = Ae q ((-w + V2uJe)tj + Be q ((-u - V2uJe)t 
Expanding this solution in terms of e, 

x(t) _ ,ff (-")" + "(^(-")-- + -) ),.. 



»E 



(-w)" - n(V2uje{-u) n - 1 + ...) 



(-«) 



(^)E^" + ^- s )^£tSt( 



(17) 



-Lj) n -H n + ... 

(A + B)zi(t) + (B - A)W — i 2 (*) + ..., (18) 

V co 

in zero approximation we get the first solution 

Xl (t) = e q (-ojt). (19) 

In the linear approximation we obtain the second solution in the form 



x%{t) = t—e q (-U)t). 



(20) 



In order to prove that solutions x\{i) = e q (—u>t) and X2(t) = t-4e g (— wt) are 
linearly independent, we check the q- Wronskian : 



W q = 



e q (-wt) tj t e q (-ujt) 

D q (e q (-U)) D q (tf t e q (-ut)) 



, or 

W q = -we,(-«t) (e g (-U) - t—e q (-ojt)J . (21) 

Here we show that the term in parenthesis is not identically zero. For q > 1, by 
using the infinite product representation of e q {x) (??), we get 

<*(-«*)= n (i - (i - J) £«*) . ( 22 ) 



t-jrlne<,(-wt) = Y^ 



or 



where 



V 9/9" 



t—e q (-ujt) = Ae g (-ut), 



^E 



00 — u 1 — , 

9/9 



■■ nl-d-fj^Wt 



Expanding the denominator, we have 



n 1 \^ / iV 1 



^ SlH'-y^JSl 1 ^)^ 



00 / 

- M 1 - 

1=0 v 


1 \ l+1 °° 1 

y 7 n=0 y 


00 ((1 


[Z] 1 - q 


where \t\ < -2-. We know that 

1 1 w 






9 Q OO / 



/=1 
Then we rewrite 



A=-^-Ln q (l-(l-- q y, ! ). 1 2:, l 



(23) 



ln(l _ a;) = -a; - — - — -... = - V - 
V ; 2 3 ^ Z 

z=i 

and the g-analogue of this expression is given as [6] 

00 1 
Ln g (l-x) = -J2w- ^4) 



The second solution X2(t) can also be rewritten in terms of q- logarithmic 
function 

x 2 (t) = t—e q (-ut) = —jLn q (l - ( 1 - - J utj e q (-ut), (26) 

where \t\ < &-. Finally, the g-Wronskian is not vanish 

W q = -w (e q (-Lot)) 2 (l - -^-jLn q (l - M - - J ojtjj £ 0, 

since the term 

1 Ln a ( 1 - ( 1 )ut 

couldn't be identically zero. 

We can also rewrite this in terms of g-logarithm, which instead of linear in 
t term for q — 1 case, now includes infinite set of arbitrary powers of i, 



(i - !)«*) 



i ^{^- q )^ L . . 
x 2 (t) = —Y,- — \j^ e ^-^ 

= -Y—Im q (l-(l-~Jut\e q (-ut). (27) 

It is easy to check that for q — > 1 our solution reduces to the standard second 
solution ie _LJ *. 

Combining the above results we find the general solution in the degenerate 
case as 

x(t) = Ae q (-ut) + Bt—e q (-cut). (28) 

In Figure [7] and Figure [8] we plot particular solutions with constant (A = 
B = 1) and with q-Periodic modulation, respectively. 

In Figures l9l and fTUl we plot solution with g-Periodic function modulation at 
different small scales. Comparing these figures we find very close similarity, this 
why g-periodic function modulation leads to the self-similarity property of the 
solution. 



5 Degenerate Roots for Equation Degree N 

g-Damped oscillator considered in first section is an example of constant coef- 
ficient q-difference equation of degree two. The problem of radiation damping 
leads to a constant coefficient equation of degree three pQ. Here we consider 
generic constant coefficient q-difference equation of degree N. Then the result for 
degenerate roots obtained in previous section can be generalized to this equation 
of an arbitrary order. 
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Figure 7: Critical case 




Figure 8: Critical case with periodic function 
The constant coefficients ^-difference equation of order N is 



J2 a k D k x{t) = 0, 



fc=0 



5> fc A fc =0. 



fc=i 



(29) 



where au are constants (or (/-periodic functions). By substitution 

x = e q {Xt) (30) 

we get the characteristic equation 

N 



It has N roots. Suppose (Ai, A2, ..., Xn) are distinct numbers. Then, the general 
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Figure 9: Self-similar micro structure at scale 0.5 



solution of (1291) is found in the form 



N 



c(t) = y^^c k e q (\ k t). 



(31) 



fe=i 
In case, when we have n-degenerate roots 

{D + uj) n x = Q, 

by substituting (|3U)l . characteristic equation is found as 

(X + oj) n = 0. 

Then the linearly independent solutions for these degenerate roots we can obtain 
in the following form : 

xi{t) = e q (-ut), 

%2{t) = t—e q (-Lot) = — -Ln q ( 1 - (1 )ut J e q (-ut), 



x 3(t) = ( t— ) e g (-ut) 



or by using the commutation relation [i, 4A = — 1, up to linearly dependent 
solution, it can be written as 

d 2 

«n(*) = *" ^S=T e «(- W *)- 

For construction solutions with degenerate roots we need following proposi- 
tions : 
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Figure 10: Self-similar micro structure at scale 0.05 



Proposition 5.0.1 We have following commutation relation 

= -D (32) 



4e D 



implies 



d d 

t—D = D [t 1 

dt \ dt 



Proof 5.0.2 By definition of D q operator the commutation relation can be 
found as follows 



>> 



= t ±(fJMLzM). D ( t S.) 

dt V (q- 1)* / V dtj 

(q - l)t (qf'{qt) - /'(*)) -(«-!) (/(<?*) - /(*)) \ ?*S " */'(*) 



f 



"Af 



fe - l) 2 i 2 



(9-1)* 



(33) 



which implies 



'> 



= -D. 



Proposition 5.0.3 



t A D " = D n ( t ± _ n 

dt \ dt 



(34) 



Proof 5.0.4 By using mathematical induction: 

For n = 1, /rom i/ie above commutation relation it is easy to see. 



12 



Suppose it is true for n: t-^D n = D n (t^ - n) 
And we should show that it is true for n + 1 



t-D n+1 = t-D n D = D n (t--n)D = D n t-D-nD n+1 
dt dt \ dt J dt 

= D n D (tj- - 1 ) - nD n+l = D n+1 (tj - (n + 1 ! ) ( :r> i 
Proposition 5.0.5 We have more general relation in the following form 

t—(uj + D) n = {uj + D) n t^- - n(u + DY^D. (36) 

Proof 5.0.6 By using mathematical induction : 
For n = 1 , 

d , . d d d ( d 

^ + ^ = t d t U + t dt D = t dt" + D { t dt- 1 

= («+*)*!- j> 

Suppose this relation is true for n: 

tj(uj + D) n = (w + D) n tf- - n(u + D) n ' l D. 
Now we prove that it is true for n + 1 

t^-(oj + D) n+1 = t^-{u + D) n (u + D) 
dt dt 

= ({uj + D) n t^--n{uj + D) n - l D\{uj + D) 
= (u + D) n Uuj + D)tj - Dj -n(uj + D) n D 

= {u + D) n+ Hj-{n + l){u + D) n D 

Using the operator identity (|36|) we can show that if x$ is solution of 
(D + w)x Q = => (D + u) 2 x = =*• ... =► (D + u) n x = 0. 



t±(D + w) n z =0 

at 

(D + u}) n t—x - nDiuj + D) n - l x = 
dt 

{lu + D)"- l x = 0^>(lu + D) n Xl = 0, 
13 



where 
Then, 

is solution of 
and as follows 
e.t.c. And then, 

is solution of 



d 
X\ = t — x . 
at 



d 
at 



{D + ujf Xl =0, 
(D + uj) n Xl =0, 



X n -1 — ^n-1* 



(£> + w)"x„_i=0. 



It provides us with n linearly independent solutions xq, Xi, ..., x n -\ of TV-degree 
equation with n-degenerate roots 



(D + uj) n x = 0. 



6 Conclusions 



In conclusion we like to mention relation of our q-damped oscillator problem 
with nonlinear q-diffcrcncc problem. By following substitution 

equation (1121) leads to nonlinear the q-Riccati equation 

D q y{t) + y(qt)y(t) + Ty(t) +lu 2 =0, (38) 

so that every solution of the first one produces a solution of the second one. 

This means that (| 12[) gives linearization of the nonlinear q-difference equation 

(|38[) . Similar situation we encounter in the q-Burgers equation linearized by 

q-Cole-Hopf transformation in therms of the q-heat equation [5], 

Finally some comments about degenerate limit of our system. If in addition 

to our q-damped oscillator model (fT!2)) we consider standard oscillator with q- 

derivative friction 

d 2 
m-px{t) + jD q x(t) + kx{t) = 0, (39) 

then in the case of strong damping 7 >> m, or the small mass, both equations 
reduce to the first order q-diffcrcncc equation 

jD q x{t) + kx(t) = 0. (40) 

14 



Under the reduction, the two-dimensional phase space for the second order 
system turns into a one-dimensional one for the first order system. The first 
order system is called the "degenerate system" [7]. An arbitrary initial value 
problem in general does not apply to the degenerate system (|4"0")) j5]. For the 
second order system at t — we can attach an arbitrary value for coordinate 
x and related velocity x or D q x. However, we can describe the same physical 
system by (|40|) only after some time interval and, moreover, in such a case x or 
D q x cannot be arbitrary since D q x is completely determined by given coordinate 
x, according to PU|). When the mass m is going to zero, transition from a 
state incompatible with (|4"0)) to a compatible one is very fast. Acceleration at 
the initial time is very high (the related velocity is changing very fast). The 
transition to the massless limit can be well approximated by the discontinuous 
jumping condition like in [7J: the energy of the system cannot be changed by a 
jump. The jumping condition implies that under the jumping the coordinates 
of the system remain invariant and only the velocities can be changed. Details 
of this study we present in our future work. 
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